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ABSTRACT 



Context. The M16 nebula is a relatively nearby Hii region, powered by O stars from the open cluster NGC 6611, which borders to a Giant 

Molecular Cloud. Radiation from these hot stars has sculpted columns of dense obscuring material on a few arcmin scales. The interface 

between these pillars and the hot ionised medium provides a textbook example of a Photodissociation Region (PDR). 

Aims. To constrain the physical conditions of the atomic and molecular material with submillimeter spectroscopic observations. 

Methods. We used the APEX submillimeter telescope to map a ~ 3' x 3' region in the CO J = 3-2, 4-3 and 7-6 rotational lines, and a 

subregion in atomic carbon lines. We also observed C'*0(3-2) and CO(7-6) with longer integrations on five peaks found in the CO(3-2) map. 

The large scale structure of the pillars is derived from the molecular lines' emission distribution. We estimate the magnitude of the velocity 

gradient at the tips of the pillars and use LVG modelling to constrain their densities and temperatures. Excitation temperatures and carbon 

column densities are derived from the atomic carbon lines. 

Results. The atomic carbon lines are optically thin and excitation temperatures are of order 60 K to 100 K, well consistent with observations 
of other Hii region-molecular cloud interfaces. We derive somewhat lower temperatures from the CO line ratios, of order 40 K. The Ci/CO 
ratio is around 0.1 at the fingers tips. 
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1. Introduction 

The interplay of massive stars with their surrounding interstel- 
lar medium, though of critical importance in the energy budget 
of a galaxy, is still poorly understood due to observational lim- 
itations, and to the generally large distance to high mass star 
forming regions. Located only 1.8 + 0.5 kpc away (Bonatto 
et al. l2006t . the M16 (Eagle) nebula is one of the best tem- 
plates for detailed analysis of the environment of high-mass 
stars. This nebula is associated with the NGC 6611 star clus- 
ter, which contains more than two dozen stars of spectral type 
earlier than BO (Duchene et al. 200l) and hundreds of lower 
mass stars and has an age of 1.3 + 0.3 Myr (Bonatto et al. 
I2006> . To the south of this cluster, the radiation from these hot 
stars has sculpted columns of dense obscuring material on a 
few arcmin scales, usually referred to as 'fingers' or 'pillars', 
of which Hester et al. |1996 1 presented stunning Hubble Space 
Telescope (HST) images. The boundary between the Hii re- 
gion and the dense molecular gas provides an archetype of a 
Photodissociation Region (PDR). 

Indications of present-day star formation near the tips of the 
fingers are seen e.g. at infrared (IR) wavelengths (Thompson 
et al. 2002). The fingers have been mapped in the J = 1 - 



line of CO and isotopologues by Pound fl998') with the BIMA 
array. White et al. (1999) used the JCMT 15 m and OSO 20 m 
radiotelescopes to map the same region in various molecular 
lines and in the submm continuum. They report CO(2-l) and 
CO(3-2) peak brightness temperatures above 40 K and 60 K, 
respectively, and a factor 2-3 lower temperature in CO(l-O). 

Here we report on new observations in two submillimeter 
windows, conducted with the APEX' 12 m telescope. We have 
mapped an area covering the three pillars in CO rotational lines. 
In addition, we have mapped the fingertips in the ^^Pi-^Pq and 
^P2-^Pi transitions of atomic carbon to probe the PDR. Details 
on the observations are given in Sect. |2l Results of modelling 
the dense molecular phase are addressed in Sect. |3j while the 
properties of the PDR itself are discussed in Sect.|4] 



2. Observations 

We have used the dual colour FLASH heterodyne receiver 
(Heyminck et al., this volume) to simultaneously observe the 
J = 4 - 3 and J = 7 - 6 lines of CO, as well as, also simulta- 
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neously, the 492 GHz and 809 GHz transitions of atomic car- 
bon. The CO J = 3 - 2 Hne was mapped using the APEX-2a 
(Risacher et al., this issue) facility receiver With both instru- 
ments, the backend used was the MPIfR Fast Fourier Transform 
Spectrometer (Klein et al., this volume), which provides 16 384 
spectral channels over a 1 GHz bandwidth. Calibration was 
done by observing the sky, hot and cold loads and using an 
atmospheric model. This provides spectra in antenna tempera- 
ture unit. Efficiencies have been determined on observations of 
planets for all receivers (Giisten et al., this issue) and applied 
to the measured signal. The spectra were then converted to the 
main beam temperature scale, using beam efficiencies of 0.73 
at 345 GHz, 0.60 at 460 and 492 GHz, and 0.47 at 809 GHz. 
The overall calibration uncertainty of the data is estimated to 
be ~20%. The pointing was checked every hour on the nearby 
sources Sgr B2(N) and GlO.62-0.38 and was found to be sta- 
ble within 3". However, parts of the observations were done 
at elevations above 80°, where the tracking becomes somewhat 
less accurate. 

We performed position switched observations, with a fixed 
reference position to the east of the easternmost finger and at 
an offset of 120", 0" relative to our (0,0) position, which is 
0^2000 = 18''18"52^19, ^2000 = -13°48'58'.'6. All offsets used in 
the following are relative to that position. We mapped a 2' x 2' 
area with 7" spacing in both directions in CO(4-3) and CO(7- 
6), and 3' x3' with 10" spacing in CO(3-2). Based on the peaks 
found in the CO(3-2) map, five positions were also observed 
in the (3-2) transition of C"*0 (see Sect.Oand Table|2}. Due 
to moderate weather conditions for CO(7-6) in the first ob- 
serving runs, we repeated the CO(4-3) and (7-6) observations 
with longer integrations in April 2006 toward these peak posi- 
tions. A limited 40" x 20" region covering the tip of the eastern 
most finger and a 30" strip on the tip of the middle finger were 
mapped in atomic carbon lines, with 7" spacing in both cases 
(Sect.|4}. The log of the observations is reported in Tabled 



3. Molecular material 

3.1. Large scale structure 

A map of the CO(3-2) emission is shown superimposed on an 
8 yum Spitzer image in Fig.^ The 8 fim emission mostly comes 
from small dust grains and polycyclic aromatic hydrocarbons 
(PAHs) excited by the incoming radiation. It traces the edge of 
the dusty structure, while the molecular emission peaks some- 
what deeper in the columns, as expected from a transition from 
ionised to atomic to molecular medium. The integrated inten- 
sity in the CO(4-3) line is shown overlaid on the HST data from 
Hester et al. J 19961 1 in Fig. |2 Given the pointing uncertainties 
mentioned above, the structure of the CO(4-3) agrees remark- 
ably well with the columns of dust seen in the HST image. 

The 3-dimensional kinematics of the region has been dis- 
cussed e.g. by Pound (1998 1. Our data are fully consistent with 
the picture that they derived. In particular, our data confirm that 
the CO emission along the easternmost column arises from two 
distinct cores. We also observed velocity gradients (see below) 
with position angles mostly pointing to the nearby O stars. 



Table 1. Log of the observations. The Tjys column gives typical 
single side band system temperatures. 



Transition 


Frequency 


T 

sys 


Date of observation 




[GHz] 


[K] 




CO(3-2) 


345.796 


230 


17-18 Nov. 2005 


C'**0(3-2) 


329.331 


330 


18 Nov. 2005 


CO(4-3) 


461.041 


800 - 


1 28 Jul. 2005 + 








[ 1 Oct. 2005 + 


CO(7-6) 


806.652 


2500 ' 


7-8 Apr. 2006 


Cl3Pi-3pn 


492.161 


850 


1 02 Aug. 2005 + 


Cl 3P2-^Pl 


809.342 


2000 


' 04-05 Oct. 2005 



-i3°4g' 



-13°50' 




1B"19'"00" 18"1B'°55' 18"18'"50" 

Fig. 1. Contour map of the CO(3-2) emission, integrated over 
the H-20 to +30 km s"' Vlsr range, overlaid on a Spitzer 
8 yum image from the GLIMPSE survey (Benjamin et al. l2003> . 
Contours correspond to 30, 60,... 180 K km s ' in integrated 
Tmb- The APEX beam at 345 GHz is shown in the upper left 
corner. The dotted line shows the limits of our CO(3-2) data. 
The asterisks show the positions where we ran LVG modelling 
(Sect. 13. 3> and the frames drawn in dashed lines delineate the 
areas mapped in atomic carbon lines (Sect.|4}. 



3.2. Dense cores 

On the tip of the eastern most column, we measure a differ- 
ence in Vlsr of 2.15+0.07 km s"' between offset positions 
(+10",+20") and (0,-10") in the CO(3-2) line. This corre- 
sponds to a velocity gradient of 7.8+0.3 km s"' pc"', in good 
agreement Pound (I1998> . who found a value of7.1 kms ' pc ' 
when scaled to the distance of 1 .8 kpc used in the present paper. 

We derived excitation temperatures by correcting the main 
beam temperatures for the Rayleigh- Jeans approximation, as- 
suming that the CO lines are optically thick and that the 
emission is filling the main beam. At position offset (-10",0), 
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LSR Velocity (km s ') 

Fig. 2. Left hand side: Contour map of the CO(4-3) emission, integrated over 20-28 km s overlaid on an HST composite 
image ([O iii], Hq, and [S n] filters) from Hester et al. (T996^. Contour levels are shown for JTmb<^v from 40 to 165 K km s"' 
in steps of 25 K km s The circle marks the size of the APEX beam for the CO (4-3) line and is placed at the position toward 
which the spectra in the right hand side of this figure were taken. Right hand side, top to bottom: Spectra of the CO (7-6), (4-3), 
(3-2), the C"*0 (3-2) and the Ci ^?2-^?\ and ^^Pi-^Po lines taken toward the -10", offset position marked in the left part of this 
figure. 



the three transitions of CO give similar temperatures around 
54+7 K, suggesting that the molecular material is close to ther- 
mal equilibrium. 

3.3. Gas temperature and column density 

To determine the physical parameters at different positions in 
the fingers, we used the CO observations and a spherically sym- 
metric Large Velocity Gradient (LVG) statistical equilibrium 
code, with the cosmic background as only radiation field. In 
order to compare the CO(7-6) and CO(4-3) observations with 
the LVG results, we smoothed the CO(7-6) data to the spatial 
resolution of the CO(4-3). 

From the CO(3-2)/C"*0(3-2) line ratios, we estimated the 
optical depths of the CO(3-2) Unes, and derived CO column 
densities by comparing with LVG simulations. In this calcula- 
tion, we used '^O/'^O=407, as derived from the relation given 
by Wilson & Rood JI994> as a function of galactocentric dis- 
tance, with dec - 6.3 kpc (Bonatto et al. l2006> . For positon (- 
10",0"),the non-detection of C"^0(3-2) atVLSR = 23.6kms-' 
gives an upper limit to the CO column density of 10'^ cm"^. For 
the other positions, the LVG modelling gives a range of possi- 
ble column densities, as reported in Table |2j our results agree 



Table 2. Results of the CO fines modeUing. Cols. 4 and 5 
give the gas kinetic temperatures, as derived by comparing the 
observed Tmb of CO(4-3) and (7-6), respectively, with the 
LVG models. Offsets are relative to the (0,0) position given in 
Sect. 12 



Offset 


Vlsr 


N(CO) 




TcO(4-3) 


TcO(7-6) 


[","] 


[km s-'] 


[cm-2] 




[K] 


[K] 


-10,0 


23.6 


< 1 X 10" 




45+7 




-10,0 


25.4 


8 X 10" - 


1 X 10'« 


54±8 


56±9 


-h40,-60 


24.6 


8 X lO" - 


6x 10'« 


40±6 


42±6 


-30,-70 


22.4 


6x 10" - 


1 X 10" 


50±9 


57±7 


-50,-100 


21.3 


6x 10" - 


1 X 10" 


30±4 


32±6 


0,-100 


214 


2x 10" - 


6x 10'** 


51±9 


47±7 



with the average CO column density derived by Pound J1998> 
and with the C'**0 column densities of White et al. J1999> . 

Using these values for the CO column densities, we further 
constrained the gas properties by running LVG calculations for 
several densities and temperatures. We found that the CO(7- 
6) and (4-3) lines are thermalized for densities higher than 
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Table 3. Peak main beam temperatures of the [Ci] lines. The excitation temperatures and line opacities are computed after 
Zmuidzinas et al. ( 1988 1, using the Tmb scale for both lines. Typical uncertainties are of order -20 K to +30 K for Tex, and +0.1 
for opacities (see text). Column densities are derived from Tgx and t(1-0) after Keene et al. J1985> . Offsets are relative to the 
(0,0) position (Sect.E). 



Offset Tmb(I-O) vlsr Tmb(2-1) vlsr Te. t(1-0) t(2-1) N(Ci) 

["] ["] [K] [kms-'] [K] [kms-'] [K] [cm"-] 



-21 





10.8 


25.37 (0.02) 


11.6 


25.94 (0.06) 


61 


0.24 


0.31 


3.2x10'' 


-14 





12.2 


25.54 (0.03) 


14.3 


25.76 (0.09) 


71 


0.23 


0.31 


3.8x10" 


-7 





9.2 


25.60 (0.06) 


10.4 


25.70 (0.17) 


66 


0.18 


0.24 


2.9x10" 








7.2 


25.75 (0.06) 


9.7 


25.52 (0.10) 


93 


0.09 


0.14 


1.8x10" 


-35 


-63 


10.3 


23.10(0.06) 


14.3 


22.95 (0.12) 


103 


0.12 


0.19 


2.6x10" 



a few 10"* cm Both lines are optically thick at the column 
densities that we computed (Table EJ; therefore, the observed 
Tmb give direct measures of the kinetic temperature of the gas. 
Other authors (Pound 1998 White et al. 1999 l reported den- 
sities of 3 - 5 X lO'* cm""* in the fingers, with higher values 
(~ 2 x 10^ cm"^) at the fingertips. Thus, the temperatures re- 
ported in TablelJlare derived from LVG models with densities 
above lO** cm"^. For lower densities, they would correspond to 
lower limits to the true kinetic temperature of the gas. Our esti- 
mates are computed by assuming beam filling factors of 1 . The 
temperatures we derived do not change significantly over the 
range of column densities computed for each position and the 
variation is, in any case, within the error bars, corresponding to 
a calibration uncertainty of 20% for each fine. 

4. Atomic carbon lines 

In the north-south direction, [Ci] 492 GHz extends over at least 
~30", or 0.25 pc at the distance of M16 (see map in Fig. |3j- 
Using Eqs. (7) and (8) in Zmuidzinas et al. (Jj)88 1, we derived 
excitation temperatures and opacities for both [Ci] lines, us- 
ing the measured main beam temperatures for both transitions. 
This assumes a uniform excitation temperature, and that the 
levels are thermally populated. Since both lines were observed 
simultaneously, systematic errors should cancel out in this anal- 
ysis, so that using a 10% uncertainty for each measured Tmb 
is reasonable. The resulting uncertainties are of order -20 to 
+30 K for Tex and ±0. 1 for opacities. 

The results are shown in Table |3l for the central part of the 
map covering the tip of the first finger, and for offset position 
(-35",63") at the tip of the second finger, where we have good 
S/N spectra in both transitions. We find optically thin lines, 
with opacities in the range 0. 1-0.3, and excitation temperatures 
in the range 60-100 K, i.e. in the upper part of the temperature 
range computed by Zmuidzinas et al. (J_988 1 in various sources. 
We also derived excitation temperatures from LVG modelling, 
and found similar results, in the range 80-120 K. 

We computed Ci column densities using the same formal- 
ism as described in the appendix of Keene et al. {1985}. They 
are reported in the last column of Table |3l The corresponding 
Ci/CO ratios are around 0. 1 at the tips of both fingers. However, 
Ci and C'^O emissions are likely to trace distinct material, and 
additional observations are considered to constrain the varia- 
tions of Ci/CO in the region. 



5. Summary and conclusion 

We have observed the transitions of atomic carbon for the first 
time toward the M16 nebula, and derived excitation tempera- 
tures for the atomic material around 60-100 K, similar to the 
values found in other Hii region-molecular cloud interfaces. 
The kinetic temperatures derived from CO transitions seem 
somewhat lower, around 40 K in most positions. This is consis- 
tent with the picture where the atomic carbon traces the PDR 
at the edge of the dense molecular columns. The Ci to CO ratio 
is ~0.1 at the fingers tips. Additional observations and com- 
parison with PDR models are needed to further constrain the 
atomic to molecular transition. 
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